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Abstract 

In this paper a microscopic state-of-the-art approach to calculation of the nuclear matrix element 
M^^ for neutrinoless double beta decay with an account for nuclear deformation is presented in 
length and applied for ^^Ge, ^^'^Nd and ^^'^Gd. The proton-neutron quasiparticle random phase 
approximation (QRPA) with a realistic residual interaction [the Brueckner G matrix derived from 
the charge-depending Bonn (Bonn-CD) nucleon-nucleon potential] is used as the underlying nuclear 
structure model. The effects of the short range correlations and the quenching of the axial vector 
coupling constant gA are analyzed. The results show that neutrinoless double beta decay of ^^*^Nd, 
to be measured soon by the SNO-I- collaboration, provides one of the best probes of the Majorana 
neutrino mass. This confirms our preliminary conclusion in Ref. (l7| . 
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I. INTRODUCTION 



Neutrinoless double beta decay (Oz//3/3-decay) is a second order nuclear weak decay process 
with the emission of two electrons only [1-3]: {A, Z) — )■ {A, Z + 2) + 2e~. This process 
is forbidden in the standard model (SM) of electroweak interaction since it violates the 
conservation of the total lepton number. The observation of OuPP decay will immediately 
prove the neutrino to be identical to its antiparticle (a Majorana particle) as assured by 
the Schechter-Valle theorem 0|. Thereby, Ou^P decay offers the only feasible way to test 
the charge-conjugation property of the neutrinos. In addition, the existence of Oz//3/3 decay 
requires that the neutrino is a massive particle. 

The fact that the neutrinos have non-vanishing masses has firmly been established by 
neutrino oscillation experiments (see, e.g., Ref. j5|). However, the oscillation experiments 
cannot in principle measure the absolute scale of the neutrino masses. One of the possible 
ways to probe the absolute neutrino masses at the level of tens of meV is to study Oi^PP de- 
cay. Provided the corresponding Oz//3/3-decay rates are accurately measured, reliable nuclear 
matrix elements (NME) M^" will be needed to deduce the effective Majorana neutrino mass 
from the experimental half-lives of the decay. 

One of the best candidate for searching Oz//3/3 decay is ^^°Nd since it has the second highest 
endpoint, Q^p =3.37 MeV, and the largest phase space factor for the decay (about 33 times 
larger than that for ''^Ge, see e.g. The SNO-f- experiment at the Sudbury Neutrino 

Observatory will use Nd loaded scintillator to search for neutrinoless double beta decay by 
looking for a distortion in the energy spectrum of decays at the endpoint [6|. SNO+ will 
be filled with 780 tonnes of liquid scintillator. The planned loading of 0.1% of the natural 
Nd translates into 43.6 kg of the isotope ^^°Nd. It is expected to achieve the sensitivity of 
— 5 ■ 10^^ years after one year of running, with the best final value of about 3-4 times 
longer (without enrichment of the dissolved Nd). With the NME M^" = 4.74 of Ref. 0], 
obtained within the proton-neutron quasiparticle random phase approximation (QRPA) with 
neglect of deformation, already the initial phase of SNO-I- will be able to probe m^^j ~ 100 
meV, and will finally be able to achieve sensitivity of m^^ ~ 50 meV corresponding to the 
inverse hierarchy (IH) of the neutrino mass spectrum. 

However, ^^°Nd is well-known to be a rather strongly deformed nucleus. This strongly 
hinders a reliable theoretical evaluation of the corresponding Oz//3/3-decay NME (for instance, 
it does not seem feasible in the near future to reliably treat this nucleus within the large- 
scale nuclear shell model (LSSM), see, e.g., Ref. (sj). Recently, more phenomenological 
approaches like the pseudo-SU(3) model joj, the Projected Hartree-Fock-Bogoliubov (PHFB) 
approach [loj and the interacting boson model (IBM-2) [lH have been employed to calculate 
M^" for strongly deformed heavy nuclei (a comparative analysis of different approximations 
involved in the models can be found in Ref. [l2|). The results of these models generally 
reveal a substantial suppression of M^" for ^^°Nd as compared with the QRPA result of 
Ref. [tI where ^^°Nd and ^^''Sm were treated as spherical nuclei. However, the calculated 
NME M^'^ for ^^°Nd reveal rather significant spread. 

The most microscopic way up-to-date to describe the effect of nuclear deformation on 
PP-decaj NME M^" and M'^" is provided by the QRPA. In Refs. a QRPA approach 

for calculating M'^'^ in deformed nuclei has been developed. Theoretical interpretation of 
the experimental NME M^^p which have been obtained for a dozen__of nuclei [l3| provides a 



test of different theoretical methods. It was demonstrated in Refs. |l4j-|l6j that deformation 
introduces a mechanism of suppression of the M'^" matrix element which gets stronger when 
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deformations of the initial and final nuclei differ from each other. A similar dependence of 
the suppression of both M^'^ and M"^ matrix elements on the difference in deformations has 
been found in the PHFB [loj and the LSSM ^. 

In the previous rapid communication [13] we reported on the first QRPA calculation of 
jy^ou i50]\j-(^ with an account for nuclear deformation. The calculation showed a suppres- 
sion of M^'^ by about 40% as compared with our previous QRPA result for ^^ONd [7| that 
was obtained with neglect of deformation. In this paper we give the details of the calcula- 
tion of Ref . [l^ , and also include different nuclei in the analysis. In addition, the effects of 
the short range correlations and the quenching of the axial vector coupling constant qa are 
considered. Making use of the newest NME, one may conclude that 0^(3(3 decay of ^^°Nd, to 
be searched for by the SNO-I- collaboration soon, provides one of the best sensitivities to the 
Majorana neutrino mass and may approach the IH region of the neutrino mass spectrum. 



II. FORMALISM 

For the light Majorana neutrino exchange mechanism, the inverse Oz//3/3-decay lifetime is 
given by the product of three factors. 



(T,%)-^ = G°^|M'°f mj,, (1) 

where G^'^ is a calculable phase space factor, M'^" is the Ou/S/S nuclear matrix element, and 
mpjs is the (nucleus-independent) "effective Majorana neutrino mass" which, in standard 
notation [isj, reads 
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(2) 



1=1 

rrii and Uei being the neutrino masses and the mixing matrix elements, respectively. The 
NME includes both Fermi (F) and Gamow- Teller (GT) transitions, plus a small tensor (T) 
contribution j2[, 

V1.25/ Y / 

In the above expression, is the effective axial coupling in nuclear matter, not necessarily 
equal to its "bare" free-nucleon value g^ — 1.25. We note that for g^ = 1-25 nuclear matrix 
element M'^" coincides with the standard definition M^" = -M'^p"^ + M^^ + M^" . With 

the conventional prefactor oc g\ in Eq. ([3]), the phase space G'^" becomes independent of gA- 
The NME M^^ for strongly deformed, axially-symmetric, nuclei can be most conveniently 
calculated within the QRPA in the intrinsic coordinate system associated with the rotating 
nucleus. This employs the adiabatic Bohr-Mottelson approximation that is well justified for 
^^°Nd, which indeed reveals strong deformation. As for ^^"^Sm, the enhanced quadrupole 
moment of this nucleus is an indication for its static deformation. Nevertheless, the exper- 
imental level scheme of ^^°Sm does not reveal a clear ground-state rotational band. In this 
work we treat ^^°Sm in the same manner as ^^°Nd. However, a more elaborated theoretical 
treatment going beyond the simple adiabatic approximation might be needed in the future 
to describe the nuclear dynamics of this nucleus. 

The adiabatic, or the strong coupling (see, e.g., the monograph [l^), approach fails, 
however, for small deformations since the Coriolis force gets large and mixes states with 
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different K. Nuclei ^^Ge and ^^Se have rather small deformations and the so-called weak 



coupling, or no alignment, limit [19| seems to be more suitable. In this limit the Coriolis 
force becomes so strong that the angular momenta of the valence nucleons get completely 
decoupled from the orientation of the core. Such a case would deserve a separate detailed 
study, and the adiabatic approach to description of excited states of deformed nuclei is 
adopted in the present application of the QRPA with a realistic residual interaction. Nev- 
ertheless, one might already anticipate without calculations that in the weak coupling limit 
the calculated observables should reveal smaller deviations from the ones obtained in the 
spherical limit than those calculated in the strong coupling limit of the present work. In 
this connection it is worth to note that spherical QRPA results can exactly be reproduced in 
the present calculation by letting deformation vanish, in spite of the formal inapplicability 
of the strong-coupling ansatz for the wave function in this limit. 

Nuclear excitations in the intrinsic system \K'") are characterized by the projection of 
the total angular momentum onto the nuclear symmetry axis K (the only projection which 
is conserved in strongly deformed nuclei) and the parity vr. In Ref. [l6| the structure of the 
intermediate lO"*") and ll"*") states was obtained within the QRPA to calculate 2z//3/3-decay 
NME M'^^ . Here, the approach of Ref. [l^ is straightforwardly extended to calculate all 
possible \K'^) states needed to construct the NME M^" . 

The intrinsic states \K'^^m) are generated within the QRPA by a phonon creation oper- 
ator acting on the ground-state wave function: 

\K'" ,ra) = <5m,_ft'|0g^s.) j Qln,K = -^pn,K^pn,K ~ ^pn,K^pn,K ■ (4) 

pn 

Here, A^^ ^^ = a],aj^ = ctpCtn are the two-quasiparticle creation and annihilation 

operators, respectively, with the bar denoting the time-reversal operation. The quasiparticle 
pairs pn are defined by the selection rules Qp — Qn = K and Tip^n = vr, where vr,- is the single- 
particle (s.p.) parity and Qr is the projection of the total s.p. angular momentum on the 
nuclear symmetry axis (r = p,n). The s.p. states \p) and |?7.) of protons and neutrons are 
calculated by solving the Schrodinger equation with the deformed axially symmetric Woods- 



Saxon potential [16[. In the cylindrical coordinates the deformed Woods-Saxon s.p. wave 
functions IrQr) with Qr > are decomposed over the deformed harmonic oscillator s.p. 
wave functions (with the principal quantum numbers {NuzA)) and the spin wave functions 
|E = ±i): 

|rfi.) = J2 bNn.^\NnzA^ = n,-J:)\J:), (5) 

where N = n± + {n± = 2np + |A|), and are the number of nodes of the basis 
functions in the z- and p-directions, respectively; A = — S and S are the projections 
of the orbital and spin angular momentum onto the symmetry axis z. For the s.p. states 
with the negative projection Vl^^ = that are degenerate in energy with Q^- = \ flr\, the 



time- reversed version of Eq. ([5]) is used as a definition (see also Ref. [16|). The states (t, f ) 
comprise the whole single-particle model space. 

The deformed harmonic oscillator wave functions \NnzA) can be further decomposed 
over the spherical harmonic oscillator ones |n,,/A) by calculating the corresponding spatial 
overlap integrals A'^j^^^^ = {nrlA\NnzA) (n,. is the radial quantum number, / and A are the 
orbital angular momentum and its projection onto 2;-axes, respectively), see Appendix of 
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Ref. 16| for more details. Thereby, the wave function ([5]) can be reexpressed as 

\rn,) = J2B;,\v^r), (6) 

V 

where IrjQr) = Z^C'/o^ v i vl'^r-^^ = Qr — is the spherical harmonic oscillator wave 

function in the j-coupled scheme [r] = [rirlj)), and = J^Cf^ _^ i ^ ^7Vn^n^-s ^a^^^s, 

with C-! ^ 1 being the Clebsch-Gordan coefficient. 
The QRPA equations: 



AiK) B{K) \ (X^\ _ /Xf . , 

-B{K) -A{K)) {y^ J ^K,m^Y^J^ 

with realistic residual interaction are solved to get the forward Xl^, backward amplitudes 
and the excitation energies cj^' and 10]^^ of the m-th K'^ state in the intermediate nucleus. 
The matrix A and B are defined by 

Apn^p'n'iK) = 6pn,p'n'{Ep + En) + gpp{UpUnUptUn' + V pV nV pi V n')Vpnp, 

Qphi^UpVnUp'Vfi' ~\~ 'VpUnVp'Ufi'^Vpn'p'n 
~gpp{UpUnVpiVn' + VpVnUp>Un')Vpnp,n> 

-QphiUpVnVp'Vn' + VpUnUp>Vn')Vpn'p'n (8) 

where Bp + Bn are the two-quasiparticle excitation energies, V^n.p'n' and Vp^j p/^/ are the p — h 
and p — p matrix elements of the residual nucleon-nucleon interaction V, respectively, Ur 
and Vr are the coefficients of the Bogoliubov transformation. 

As a residual two-body interaction we use the nuclear Brueckner G matrix, that is a 
solution of the Bethe-Goldstone equation, derived from the Bonn-CD one boson exchange 
potential, as used also in the spherical calculations of Ref. 0]. The G matrix elements are 
originally calculated with respect to a spherical harmonic oscillator s.p. basis. By using the 
decomposition of the deformed s.p. wave function in Eq. (jS]), the two-body deformed wave 
function can be represented as: 

b^) = E Fp,'',nJVpVn,JK), (9) 

Vp'TlnJ 

where \7]p7]n,JK) = Em^m„Cl^^j^„^Jr]pmp)\r]nmn), and i^pifnr,„ 
BP^^B'^^{-l)^"~^"CjJ^^j^_^^ is defined for the sake of simplicity ((-1)-?"-^" is the phase 
arising from the time-reversed states |n)). The particle-particle Vpn^ pi^i and particle-hole 
Vpni, p'n interaction matrix elements in the representation (|8]) for the QRPA matrices A, B 
([7]) in the deformed Woods-Saxon single-particle basis can then be given in terms of the 
spherical G matrix elements as follows: 

ypn,p'n' =-2^^ ^ Bp^^nri„Fp\,n'ri„,Gir]pVnVp'Vn', J), (lO) 

J VpVnVp'Vn' 

JX' JK' 

Vpn',p'n = '^T.T. PpVp^V,.'Fp\'mn<^iVpVn'Vp'Vn,J), (H) 

J VpVnVp'Vn' 

where K' , = Qp + fin' = ^p' + ^n- 



5 



The matrix element is given within the QRPA in the intrinsic system by a sum of 
the partial amplitudes of transitions via all the intermediate states K'": 

= M^^^l^") ; M^''{K^) = (12) 



Here, we use the notation of Appendix B in Ref. [20[, a stands for the set of four single- 
particle indices {p,p' ,n,n'}, and Oa{K'^) is a two-nucleon transition amplitude via the K'^ 
states in the intrinsic frame: 

0^{K^)= Y {0^\clcn\K^mf){K''mf\K^m,){K''m,\cl,Cr.>\Oi). (13) 

mi,mf 

The two sets of intermediate nuclear states generated from the initial and final g.s. (la- 
beled by rrii and rrif, respectively) do not come out identical within the QRPA. A standard 
way to tackle this problem is to introduce in Eq. f[T^ the overlap factor of these states 
{K'^rUflK'^mi), whose representation is given below, Eq. f ITB]) . Two-body matrix elements 
s'^^^^ of the neutrino potential in Eq. (fT2|) in a deformed Woods-Saxon single-particle basis 
are decomposed over the the spherical harmonic oscillator ones according to Eqs. (iQfTTj) : 



Adef) ^ V- pJK jpJK (sph) ( J^ /..N 



VpVp' 

VnVn' 



where JT" = y/2J' + 1, and Oe{l,2) is the neutrino potential as a function of coordinates of 
two particles, with i labeling its Fermi (F), Gamow- Teller (GT) and Tensor (T) parts. 

The particle-hole transition amplitudes in Eq. ( !T3|) can be represented in terms of the 
QRPA forward X^^ and backward Y^]^ amplitudes along with the coefficients of the Bogoli- 
ubov transformation Ur and Vr 161 : 



l^K Tni\c}pCn\Q~i) = UpVnXp^^j^TT -\- VpUnYpj^''^j^-n . 

The overlap factor in Eq. (fT3|) can be written as: 

{K-mj\K-m,) = ^Kk^Kk^ ' Y;;^.Y,^^.]ni^i^{BCSf\BCS,) (16) 



liif 



Representations for TZi^i^ and the overlap factor {BCSf\BCSi) between the initial and final 
BCS vacua are given in Ref. 3|- 



III. RESULTS AND ANALYSIS 

We have computed the NME M^" for the Ouf3(3 decays ^^Ge^^'^Se, i5°Nd^^^°Sm, 

160QJ^160j3y. 
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The single-particle Schrodinger equation with the Hamiltonian of a deformed Woods- 
Saxon mean field is solved on the basis of a axially-deformed harmonic oscillator. The 
parametrization of the mean field is adopted from the spherical calculations of Refs. 0, 20, 



211 ] ■ We use here the single-particle deformed basis corresponding in the spherical limit 
to full {4:-6)huj shells. Decomposition of the deformed single-particle wave functions is 
performed over the spherical harmonic oscillator states within the seven major shells. Only 
quadrupole deformation is taken into account in the calculation. The geometrical quadrupole 
deformation parameter P2 of the deformed Woods-Saxon mean field is obtained by fitting the 
experimental deformation parameter (3 = ^/f where is the charge rms radius and Qp 
is the empirical intrinsic quadrupole moment. The latter can be derived from the laboratory 
quadrupole moments measured by the Coulomb excitation reorientation technique, or from 
the corresponding B(E2) values 22[. We take in this work experimental values extracted 
from the B(E2) values as being more accurate. The fitted values of the parameter /32 of 
the deformed Woods-Saxon mean field, which allow us to reproduce the experimental /3, are 
listed in Table [H We label these sets of parameters as "1". The spherical hmit, i.e. ^2 = 0, 
is considered as well (labeled as "0"), to compare with the earlier results of Ref. |7|. The 
adopted here procedure of fitting is more consistent than the approximate ansatz (32 = (3 
used in Ref. 

As in Refs. 0, 16, 17, 20, 21], the nuclear Brueckner G matrix, obtained by a solution 
of the Bethe-Goldstone equation with the Bonn-CD one boson exchange nucleon-nucleon 
potential, is used as a residual two-body interaction. First, the BCS equations are solved to 
obtain the Bogoliubov coefficients, gap parameter and the chemical potentials. To solve the 
QRPA equations, one has to fix the particle-hole Qph and particle-particle Qpp renormalization 



factors of the residual interaction, Eqs. [HI As in Ref. [16|, llTJ , we determine a value of Qph 
by fitting the experimental position of the Gamow- Teller giant resonance (GTR) in the 
intermediate nucleus. Since there is no experimental information on the GTR energy for 
^^°Nd, we use for this nucleus the same gph = 0.90 as fitted for ^^Ge (this value is slightly 
different from the fitted Qph = 1.15 of Ref. [l6| because of a different parametrization of the 
mean field used here). The parameter Qpp can be determined by fitting the experimental 
value of the 2z//3/3-decay NME M^^ = 0.07 MeV"^ (l3|. To account for the quenching of the 
axial- vector coupling constant qa, we choose in the calculation, along with the bare value 
qa = 1.25, also the quenched value g'^'^ = 0.75- qa = 0.94, where the quenching factor of 0.75 
comes from a recent experimental measurement of GT strength distribution in ^^°Nd (23| . 
The two sets of the fitted values of gpp corresponding to the cases without or with quenching 
of gA are listed Table H] as cases (I) and (II), respectively. Note, that the more realistic 
procedure of fittin g adopted here also give us more realistic gpp ~ 1 values as compared 
with those of Ref. [l6| . 

Having solved the QRPA equations, the two-nucleon transition amplitudes f lT3|) are calcu- 
lated and, by combining them with the two-body matrix elements of the neutrino potential, 
the total OuPP NME M^'^ f|T2|) is formed. The present computation is rather time consum- 
ing since numerous programming loops are needed to calculate the decompositions of the 
two-body matrix elements in the deformed basis over the spherical ones. Therefore, to speed 
up the calculations the mean energy of 7 MeV of the intermediate states is used in the neu- 



trino propagator. Following Refs. [7|, |20|, |21|, we have taken into account the effects of the 
finite nucleon size and higher-order weak currents are included. Recently, it was shown [2l| 
that a modern self-consistent treatment of the two-nucleon short-range correlations (s.r.c.) 
leads to a change in the NME M^'^ only by a few percents, much less than the traditional 
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TABLE I: The values of the deformation parameter of Woods-Saxon mean field for initial 
(final) nuclei fitted in the calculation to reproduce the experimental quadrupole moment (labeled 
as "1"). The spherical limit is labeled as "0". Also the fitted values of the particle-particle strength 
parameter Qpp are listed (for both cases without (I) and with quenching (II) oi qa)- The particle- 
hole strength parameter is Qph = 0.90. The BCS overlap factor {BCSf\BCSi) (fTBl) between the 
initial and final BCS vacua is given in the last column. 



initial (final) 


P2 


9pp (I) 


9pp (11) 


{BCSi\BSCf) 


nucleus 










^^Ge (76se) 


0.10 (0.16) "1" 


0.71 


0.66 


0.74 




0.0 (0.0) "0" 


0.68 


0.63 


0.81 


isoNd (isosm) 


0.240 (0.153) "1" 


1.05 


1.00 


0.52 




0.0 (0.0) "0" 


1.01 


0.99 


0.85 


leoQd (leoDy) 


0.303 (0.292) "1" 


1.00'' 


1.00 


0.74 



°As there is no experimental value of M'^'^ for ^^"Gd, we do not renormalize the p-p interaction and use 
Opp — 1- 



Jastrow-type representation of the s.r.c. does. A very similar effect is found in the present 
calculation (see below). 

We start our discussion of the calculated Oi//3/3-decay NME by a comparison of the matrix 
elements of this work obtained in the spherical limit with the previous ones of Refs. 



that provides an important cross-check of the present calculation. Though formally the adia- 
batic Bohr-Mottelson approximation is not applicable in the limit of vanishing deformation, 
it is easy to see that the basic equations (fT2|) - (fT6|) do have the correct spherical limit. 

According to Eq. (fT2l) . the total calculated Oi//3/3-decay NME is formed by the sum of 
the all partial contributions M^'^{K'^) of different intrinsic intermediate states K'^ , with 
M^''{\K\'^) = M^'^{—\K\'^). In the spherical limit the total angular momentum J becomes a 
good quantum number, and the intermediate states corresponding to a given J'^ state 
become degenerate. In addition, in this limit each projection K of an intermediate J'" 
state contributes equally to the calculated 0^(3/3- decay NME, as a consequence of rotational 
symmetry. To represent standard spherical results in the terms of the present paper, one 
has to use the following expression for the spherical partial contribution of a projection K: 
M°''{K'^) = J2 M^''{.r)/{2J + 1) (it is easy to see that having summed over all K'^ one 

J>\K\ 

obtains the total NME M^'^ = J2 M^'^{J'^)). From this representation it can generally be 

J 

expected that the smaller is \K\, the larger the corresponding partial contribution M^^{K'^) 
should be (since simply more J's contribute, and their contributions are of the same sign in 
most cases, see Q). This behavior is in fact revealed by most of the calculation results (see 
below) . 

To test our new numerical code calculating Oz//3/3-decay NME for deformed nuclei, we have 
taken the spherical limit and used in it the _spherical harmonic oscillator wave functions as 
usually done in the QRPA calculations 0,120, 2l[. The NME calculated by different codes 
are found to be in an excellent agreement. 

For the further discussion we define the following contributions to the total Oz//3/3-decay 
NME: Mp:'{l/r) and Mg.^^(l/r) are calculated by taking into account only the Coulomb- 
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like radial dependence of the neutrino potential. The total corrections AMp:' and AMq'j, 
to Mf{l/r) and Mg^(l/r), respectively, come from the effects of the finite nucleon size 
(FNS), the closure energy, and higher order weak currents. Finally, the corrections diMp' 
and SiMo'rp, respectively, come from the effects of the s.r.c. (Jastow-like s.r.c. is denoted by 
the subscript i = I, and the self-consistent Bonn-CD s.r.c. — by the subscript i = 2). Thus, 
the final total Oz//3/3-decay NME is given by 

= M°^(l/r) + AM^" + 5,M°^ (17) 

The effect of deformation and different choices of the single-particle wave functions on the 
partial contributions M^^{K'") of different K'" intermediate states to M^"{l/r) for ^^°Nd is 
illustrated in Fig. [TJ The BCS overlap factor is neglected here for simplicity. The Fermi and 
the GT contributions are shown in the left and right panels of the figure, respectively. The 
left and middle bars for each K'^ represent the results obtained for zero deformation with 
the spherical harmonic oscillator wave functions and with the Woods-Saxon wave functions, 
respectively, and the right bar represents the result calculated with the deformed Woods- 
Saxon wave functions for the finite deformations from Table [T] (case "1"). Each bar for 
K > Q represents a sum of the equal contributions of the positive and negative projections, 
M^^iK"") + M^^i-K"") = 2M^''\k''). Because the partial contributions Mf{K^) to the 
Fermi NME in the spherical limit are non-zero only for the intermediate states of natural 
parity vr = (-1)-^, the equality M^''{K^,7r = (-1)^) = M^^^i^K - = (-1)^) {K > 

0) should hold in this limit (this simply reflects equality of the contributions of different 
projections K for a given J). This equality is nicely fulfilled in our calculations, as illustrated 
by Fig. [U Also, one can see in Fig. [1] that even for rather large deformations the partial 
contributions M^'^{K'^) (calculated with neglect of the BCS overlap) do not differ much from 
the corresponding ones of the spherical calculations. The corresponding results for the other 
nuclei show very similar pattern. 

The final results for the NME for Oi//3/3 decays ^^Ge^^^Se, i^ONd^^^^Sm, ^^^Gd.^^^^T)y 
are listed in Table HI] (now the values of the BCS overlap from the last column of Tableware 
taken into account). In columns 4 and 9 the leading contributions Mpil/r) and Mgt{^/^) 
are shown. In columns 5 and 10 the total corrections AMp and AMqt to Mp{l/r) and 
Mgt(1/''") are listed. One sees that these corrections reduce the results of the leading order 
by about 20% for the Fermi part and 40% for the GT part, that agrees with the previous 



spherical QRPA calculations [2J]. In columns 6,7 and 11,12 the corrections 5i^2Mp and 
'5i,2^GT; respectively, coming from different choices of the s.r.c, are listed. In columns 8 
and 13 both the F and GT parts of the total NME M^^^ = Mf''{l/r) + AM^^" + ^sM^^^ are 
shown (we prefer here the final value of M^^ corresponding to the modern self-consistent 
treatment of the s.r.c. |2l|). Finally, in columns 14 and 15 the Oi//3/3-decay NME M'^" 
and corresponding decay half-lives (assuming m^^=50 meV) are listed. The corresponding 
/^'^-decompositions of M*^^ are shown in Figs. [2ll3l 

By inspecting Table [III and Figs. [2][3| one sees that the difference between the spherical 
and deformed results mainly come from the BSC overlap between the ground states of the 
initial and final nuclei. As for the Qpp dependence of the Oz//3/3-decay NME, it is much less 
pronounced as compared with the dependence of the amplitude of 2z//3/3 decay. A marked 
reduction of the total M'^^ for the quenched value of qa can be traced back to a smaller 
prefactor (^^/1.25)2 in the definition of M'^"" 

The strongest effect of deformation on M"'" (the suppression by about 40% as compared 
to our previous QRPA result obtained with neglect of deformation) is found in the case 
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FIG. 1: The partial contributions M^'^{K'^) of different intermediate states to M^'^(l/r) for 
I50]sj^_j.i50gj^ in the cases of vanishing and reahstic deformations. For simphcity, the BCS overlap 
factor is omitted in these results. The Fermi Mp^{l/r) and the GT MQ^{l/r) contributions are 
shown in the left and right panels, respectively. The three bars represent (from left to right) the 
results obtained with the spherical harmonic oscillator wave functions, with the Woods-Saxon wave 
functions in the spherical limit, and with the deformed Woods-Saxon wave functions for realistic 
deformations from Table HI 



of ^^"Nd. This suppression can be traced back to a rather large difference in deformations 
of the ground states of ^^°Nd and ^^°Sm. Such an effect has been observed also by other 
authors @,ll0|,[il|. The modern self-consistent treatment of the s.r.c. l2ll makes the resulting 
M'°^ = 3.34 (without quenching), even a bit larger than the NME M'^" = 3.16 of Ref. |lj 
where the influence of the s.r.c. was completely neglected. This translates to the half-life 
= 4.1 ■ 10^^ y for the effective Majorana neutrino mass {jn^p) = 50 meV (cf. with 
= 4-60 ■ 10^5 y of Ref. (Tj). In the case of quenched c/a, the half-life Tf/a = 7-1 ■ 10^^ y 
is about twice longer as a consequence of a smaller NME M'^'^ = 2.55. 



IV. CONCLUSIONS 



In this paper a microscopic state-of-the-art approach to calculation of the nuclear matrix 
element M^'^ for neutrinoless double beta decay with an account for nuclear deformation is 
presented in length and applied to calculate M^" for ™Ge, ^^°Nd and ^^"Gd. The QRPA with 
a realistic residual interaction (the Brueckner G matrix derived from the Bonn-CD nucleon- 
nucleon potential) is used as the underlying nuclear structure model. The effects of the short 
range correlations and the quenching of the axial vector coupling constant qa are analyzed 
and found to be in accord with the spherical QRPA calculations. The strongest effect of 
deformation on M^'^ (the suppression by about 40% as compared to our previous QRPA 
result obtained with neglect of deformation) is found in the case of ^^''Nd. This suppression 
can be traced back to a rather large difference in deformations of the ground states of ^^°Nd 
and ^^°Sm, that agrees with results by other authors. The preliminary conclusion of Ref. [l3| 
that neutrinoless double beta decay of ^^°Nd provides one of the best probes of the Majorana 
neutrino mass is confirmed. 
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TABLE II: Different contributions to the total calculated NME M'°'' for decays '^'^Ge-;>^^Se, 

^^°Nd— )'^^°Sm, ^^°Gd-)'^^°Dy. The BCS overlap is taken into account. In columns 4 and 9 the 
leading contributions Mp'(l/r) and McTi^/f) are shown. In columns 5 and 10 the total corrections 
AMp and AMqt to MF{l/r) and McTi^/f) are listed. In columns 6,7 and 11,12 the corrections 
SiMp and diMcT, respectively, coming from different choices of the s.r.c, are listed. In columns 8 
and 13 both the F and GT parts of the total NME ()17p are shown (we prefer here the final value of 
j^Ou corresponding to the modern self-consistent treatment of the s.r.c. 2l||). Finally, in columns 
14 and 15 the Oz^/3/3-decay NME M'^'^ 1^ and corresponding decay half-lives (assuming m^^=50 
meV) are listed. 



A 


Def. 


9A 


M(l/r) AM 6iM 52M tot 


M{l/r) AM diM 62M tot 




-^1/2' y 

(m^^=50 meV) 


76 


"0" 


1.25 
0.94 
1.25 
0.94 


-2.83 0.69 0.40 -0.08 -2.22 
-2.98 0.73 0.40 -0.18 -2.44 
-3.15 0.78 0.45 -0.09 -2.47 
-3.31 0.82 0.46 -0.10 -2.59 


5.59 -2.49 -0.88 0.18 3.27 
7.69 -3.65 -1.47 0.27 4.31 
6.37 -2.85 -1.01 0.20 3.72 
7.28 -3.15 -1.04 0.21 4.35 


4.69 
4.00 
5.30 
4.10 


7.15 
9.83 
5.60 
9.36 


150 


"0" 


1.25 
0.94 
1.25 
0.94 


-2.09 0.51 0.33 -0.06 -1.64 
-2.16 0.52 0.33 -0.06 -1.70 
-4.07 0.99 0.67 -0.13 -3.21 
-4.12 1.00 0.68 -0.13 -3.25 


4.01 -1.86 -0.72 0.14 2.29 
4.44 -2.00 -0.73 0.14 2.58 
7.35 -3.54 -1.46 0.26 4.07 
7.69 -3.65 -1.47 0.27 4.31 


3.34 
2.55 
6.12 
4.52 


0.41 
0.71 
0.12 
0.23 


160 




1.25 


-2.14 0.51 0.32 -0.07 -1.69 


4.57 -2.04 -0.71 0.14 2.67 


3.76 


2.26 
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FIG. 2: The partial contributions M^'^{K'^) of different intrinsic intermediate states to the total 
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s.r.c, is represented by the filled bar ("M(l/r) + AM + 5iM"). The upper and lower panels show 
the results corresponding to the unquenched = 1-25 and quenched qa = 0.94 (fitted values (I) 
and (II) of respectively, see Table H]). 
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